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Abstract 

A formalism is developed to enable the construction of the effective action 
and related quantities in QED for the case of time-varying background elec- 
tric fields. Some examples are studied and evidence is sought for a possible 
transition to a phase in which chiral symmetry is spontaneously broken. 
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The proper-time formalism was used a long time ago by Schwinger Q to 
compute the Green function for an electron propagating in the presence of a back- 
ground electromagnetic field. Although the formalism is general, explicit evaluation 
of the propagator, and of the associated effective action, was possible only for the 
case of fields uniform in space and constant in time. 

Over the intervening decades, attempts have been made |^] to compute the 
corrections to Schwinger's results for the case of varying fields. These take the form 



of a derivative expansion in the fields, but even the first non-trivial corrections turn 
out to be quite unwieldy, and are, moreover, restricted to fields that do not vary 
too rapidly (else the higher terms in the expansion must be included). 

More recently, there has been cause for a new look at this problem. The 
motivation is the strange results of the GSI heavy-ion scattering experiments Q, in 
which mysterious narrow peaks are seen in the energy spectra of emitted e+e- pairs. 
Among the many theoretical ideas that have been advanced, I wish to concentrate 
on one proposed explanation that the heavy ions create a very strong and 

rapidly varying electromagnetic field, which then induces a phase transition in QED 
to a vacuum in which chiral symmetry is spontaneously broken. The observed e+e- 
peaks are due to the decay of positronium-like states in the new phase of QED. 

To study this possibility, we employ a proper-time representation for the 
vaccum expectation value of -ipip, which is an order parameter for this transition. 
This representation is |||] 



oo 







(V'V) = m I dte-"" ^U{x,t) (1) 

where 



U{x,t) = tr(^x\e^^ (2) 

Here r is the proper-time (continued to imaginary values) and m is the electron mass. 
U is the trace of a quantum-mechanical matrix element for which the Hamiltonian 
is (7 • tt)^, with 
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-n-f, = p^M - eA^{x). (3) 

The dynamical degrees of freedom are the four coordinates Xf^{T), and are the 
associated canonical momenta. A^{x) is the potential that encodes the background 
field. We are working in Euclidean space, so the 7^ are Hermitian and H is positive. 
The trace in eq. (2) is over the indices carried by the 7 matrices. For later reference, 
note that H possesses a quantum-mechanical supersymmetry [||, generated by the 
charges Q± = |(1±75)(7 • tt) with 

{Q+,Q.} = H. (4) 

The proper-time expression for (ip^) incorporates correctly all the effects of 
the background field, but completely neglects the role of the dynamical photons. In 
using eq. (1), one hopes that these photons will not affect the presence or absence 
of a phase transition induced by the background field. 

The signal for the spontaneous breakdown of chiral symmetry is that the 
limit m ^ of (ipip) should not vanish. Because of the explicit factor of m in eq. 
(1), one requires the integral to diverge as m ^ 0. In fact, one easily sees |^] that 
if the large-r behavior of U{x,t) is t~2, (ipTp) will remain finite and non-zero. If 
the falloff is more rapid, (V'V') '^iH vanish, indicating that there is no spontaneous 
chiral symmetry breaking. 

For a free fermion, U ~ l/^"^ so that {'li^ip) — > as expected. For constant 
F^i,, one finds from the analytic continuation of Schwinger's results that U{x,t) is 
a function of the two Euclidean invariants F = \{E'^ + B"^) and G = E ■ B . If 
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G = 0, F ^0, U 1/t and (V-^) ^0. If G / 0, then U ~ const., which indicates 
that ("ipip) — > oo as m =— > 0. This behavior is not to be interpreted as spontaneous 
chiral symmetry breaking, however, since there is an anomaly when G ^ that 
expUcitly breaks chiral symmetry. In what follows, we shall look at cases where 
ii^ 7^ 0, i? = 0, so the anomaly is absent. 

In a recent paper, Caldi and Vafaeisefat Q have computed U{x,t) numer- 
ically using Monte Carlo simulation techniques. For this purpose, it is convenient 
to recast U{x,t) as a path integral: 

U{x,t) = tr J DXf,Te-^, (5) 
S= r dt'L{x,x), (6) 



L{x,x) = jXf^x^ + iex^A^{x) — ^a^i,F^u{x). Note the following peculiarities: 
(i) L is complex (this is a consequence of having continued to imaginary proper- 
time); and (ii) L is matrix valued. The symbol T in eq. (5) denotes r-ordering. 

Caldi and Vafaeisefat look initially at background electric field configurations 
pointing in one direction only, for which the magnitude varies in time and in the 
one spatial variable. In particular, they consider 

^=(/(x,t),0,0) (7) 

with 
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f{x,t) = eE[cosh'^{x/Ws)]~^exp{-t^/2W^). (8) 

They find, for Wg = = 3 (in units where eE = 1) that at suitably chosen 

values for x^,U{x,t) exhibits the desired t~^/^ fallofF. Although the computations 
are complicated, their method gains credence from the following observations: they 
obtain agreement with Schwinger's analytic results for the case E = const., and, 
when Wg and Wt are taken much larger or much smaller than the above values, 
the chiral symmetry breaking goes away. This is reasonable because, for large Wg 
and Wt the configuration approximates a constant field, for which chiral symmetry 
is unbroken, whereas for small Wg and Wt the field is varying so rapidly that the 
vacuum does not have time to realign (i.e., the "sudden approximation" is valid). 
In later work Q, Caldi and Vafaeisefat have studied more realistic configurations 
involving all the spatial variables, and they continue to see chiral symmetry breaking 
for suitable values of the parameters. 

Even assuming the utter reliability of these results, one is still left with 
virtually no intuition or insight concerning the mechanism whereby chiral symmetry 
is broken. It is therefore of interest to explore these questions in a more analytic 
fashion. To this end, we look at a configuration even simpler than that chosen by 
Caldi and Vafaeisefat, to wit one in which the electric field depends only on time 
(which is, of course, Euclidean time, and which we call xq): 

E = {f{xo),0,0). (9) 
We take the associated vector potential to be 



5 



Aq = -xif{xo) 



(10) 



and 



1=0. 



(11) 



The Lagrangian then reduces to: 



L = 



-(xq + xj) + iexoXif{xo) + e<Toi/(xo) + -:{xl + x|). 



(12) 



Note the following: 

(a) For proper normalization (i.e. to obtain the known result when /(xq) = 0) we 
must take 



(b) The X2 and X3 integrals are then trivial, yielding a factor of ; 

(c) The xi integral is almost trivial, since it can be reduced to a Gaussian by a shift 
in xi; 

(d) the T-ordering is superfluous because there is only one non-trivial matrix, aoi. 
Furthermore, the trace is reduced to a summation over the eigenvalues of croi, i.e. 
to the operation 2J2a=±i (where we have abbreviated aoi by a). 

(e) The path integral must be evaluated subject to the boundary condition Xju(O) = 
Xju(r) = X/j,, where x^ is the argument of ipip. 




(13) 
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Making use of standard manipulations, one obtains 



Here F is defined by ^ = f{xo), and for any function $(r), we have defined 
($) = ^ Jq dT/^{T'). Note that U is invariant under F{xo) F{xo) + C, as it must 
be. Also note that the (F)^ term is a non-local interaction. 

We can re-express U{x,t) in terms of a local action at the cost of introducing 
an ordinary integral over a parameter A. Using 



- r dAe-"^'+^^ (15) 

TT J-oo 



with a = T and /? = 2er (F) , we have 



Uix,T) = E /^^ /l^^^O^"-^" ^^^^ 



CT=±1 



where 



Thus we are summing and integrating over a family of one-dimensional quantum- 
mechanical models defined by the Hamiltonians 



H^=p^ + V^ixo), (18) 
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V^ = [eF{xo)-\f + eaf{xo). (19) 

Here we see that H± have the standard form Q characteristic of quantum-mechanical 
supersymmetry, 

H±=p'^ + W'^ ± dW/dxo, (20) 

with 

W = eF-X. (21) 
For the purpose of quantitative analysis, it is convenient to re-express U {x, r) 

as 



1 

U{x, t) = J2 d\<x \ e-"^^ I X >, (22) 

47r r ^^^^ J-oo 

and then to insert this in eq. (1), and perform the r integral after division by m 
and differentiation with respect to m?. One obtains 



^^^^ ^ d ^ < il^jjjx) > j ^ _}_y~ f ^x<x\ ^ I X > (23) 

Therefore, we wish to compute the Green function Ga{x,x') =< x \ \ x' >, 

which obeys the equation 



8"^ 

[-^ + V^{x) + m^]G{x,x') = 5{x - x'). (24) 
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A standard expression for G is 



Gix,x') = 'ilj^{x)^^{x')e{x - x') + ^^ix')'ip^{x)e{x' - x) (25) 
where each ip obeys the homogeneous equation 

[-^ + ^- + "^']V' = 0, (26) 

subject to the boundary condition that ijjy {ip<) is weU-behaved as x ^ oo(x — oo), 
and where the Wronskian condition 



is imposed. We then have 



1 

^^""^ = 4^ y-oo"^^ ^ CWV'<(^). (28) 

Our computational strategy is to choose f(x) so that tp^ and ^< can be computed 
explicitly [p!o| , [Tl| , to insert them in eq. (28), and to determine therefrom the 
behavior of I{m) as m — > 0. If (V'V') indeed tends to a finite, non-zero value, we 
expect /(m) ~ 1/m^. Any less singular behavior will be evidence that (ipip) is 
tending to zero. 

Actually, without any further computation, we can infer from eq. (23) that 
I{m) will probably behave as 1/m?, provided there is some range of A for which 
supersymmetry is unbroken. Under these circumstances, one of H± will have zero 
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as its lowest eigenvalue, and hence the matrix element expanded in energy eigen- 
states will have a term that goes as l/m?. For the range of A (if any) for which 
supersymmetry is broken, the situation is worse: Both Hj^ and will have strictly 
positive ground state energies, so the contribution to I{rn) will be non-singular. It 
is hard to imagine a system for which the desired l/m? singularity might appear. 

As an illustrative example, we can choose 

eF{x) = jtanh(3x. (29) 

This yields a model that is exactly solvable quantum mechanically One finds 
that supersymmetry is unbroken for | A/7 |< 1, and broken when | A/7 |> 1. 
For A = ±7, the zero-energy eigenstate is not isolated but sits at the bottom of 
a continuous spectrum. Some of the energy levels of this model as a function of 
A are illustrated in Figure 1. For any values of 7 and A one can solve for ipy and 
'(/'< explicitly in terms of hyper geometric functions 2-^1- We do not reproduce the 
formulas here, since they are complicated and not particularly illuminating. When 
inserted into eq. (28), they yield the expected l/m^ behavior, i.e. no evidence for 
chiral symmetry breaking. 

It is possible to study other exactly solvable quantum mechanical models as 
well. Examples are available for which supersymmetry is unbroken for all A and 
there are others for which supersymmetry is broken for all A except A = 0. In all 
the cases we have examined, the singularity at small m of I{m) is no worse than 
1/m^. 

This result is not in conflict, of course, with the numerical work of Caldi 
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and Vafaeisefat, since their field configurations depend on at least two variables. 
In deciding how to proceed, one can think of a number of possibilities: (i) extend 
the search among the one-dimensional models in the hope that an as yet undis- 
covered class will yield the sought-for behavior; (ii) introduce new analytic 
techniques that will enable one to study the two-variable case. This will permit 
direct comparison with the Monte Carlo results; (iii) find a way to extract the small 
m behavior of I{m) (or equivalently the large t behavior of U (x, r)) without first 
having to compute the full functional forms of I{m) or U{x,t). This would lead 
to enormous simplifications not only of the analytic work but also of the Monte 
Carlo calculations, where the large r behavior is extracted by computing U (x, r) 
for several values of r and finding the slope of the best-fitting straight line on a 
log-log plot. 

As new data from GSI and Argonne are reported, one expects the relevance 
of the ideas upon which the present work is based either to wax or to wane. If the 
former, it will be interesting to see whether new insight can in fact be gained about 
the mechanism whereby time- and space-varying background fields induce a chiral 
phase transition in QED. 
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Figure Caption 

The first two energy levels of the F = jtanh/dx model discussed in the text, = 
for A < Ao = 7, and Eq = Ec = {^y — A)^ for A > 7. Ec is the energy at which the 
continuous spectrum begins. The first excited bound state exists for A < Ai = 
(7 ~ /3)^/7 (provided /3 < 7), and is given explicitly by Ei = 
/3(27 — — A^/(7 — /9)^). The diagram is symmetric for A — > —A. 

Figure available by request (chodos@yalph2.physics.yale.edu). 
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